Depinning transition for a screw dislocation in a model solid solution 
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On the basis of the classical dislocation theory, the solid solution hardening (SSH) is commonly 
ascribed to the pinning of the edge dislocations. At the atomic level, the theoretical study of the 
dislocation cores contrasts with such a prediction. Using the static molecular simulations with some 
interatomic effective potentials, we demonstrate numerically that the critical resolved shear stress 
associated with a screw dislocation in a random Ni(Al) single crystal has same order as the edge one. 
Such a result is imposed by the details of the dislocation stacking fault and the core dissociation 
into Shockley partials. The SSH statistical theory is employed to tentatively predict analytically 
the data acquired through our atomistic simulations at different Al concentration. 

PACS numbers: 62.20.F-,83.60.La 



I. INTRODUCTION 

In the alloy manufacturing the soUd sohition hardening 
(SSH) is a standard process which allows to increase the 
yield stress of a material by dispersion of some atomic- 
sized obstacles across the dislocation glide. The choice of 
the impurity and which proportion is required is an im- 
portant issue in the commercial alloy design.— On the 
condition that these impurities remain in solution, fa- 
vored either by a thermal treatment or the alloy thermo- 
dynamics, the microstructure is unchanged but the dislo- 
cations are pinned by the randomly distributed obstacles. 
The dislocation pinning yields an increase in the material 
strength without involving large inhomogencities such as 
inclusions or grain boundaries, by contrast to the other 
methods as the precipitation strengthening or the strain 
hardening. 

The dislocation impinging on a random distri- 
bution of obstacles is a standard problem of the 
theoretical metallurgy^iii^i^ii^i^iiS and the statistical 
physicsiiii^ii^ii^ii^ as well. Only recently this problem 
could have been addressed with some three-dimensional 
atomistic simulationsi^ii2ii^iiSi2ai^i2^i2^i^i2^ that shed a 
new light on points that were still a matter of debate in 
material science. However the difficulty of developing re- 
liable inter-atomic potentials for modeling dislocations in 
alloys confines the atomistic simulations to only few sys- 
tems. The Ni(Al) 7-phase is one of them and corresponds 
to the prototypical case for the binary substitutional al- 
loys with a high order energy and a high solubility limit, 
about Cai = 10 at. % which involves a broad concentra- 
tion range for the stability of the solid solution. 

The static atomistic simulations allow to compute 
the Critical Resolved Shear Stress (CRSS) for an iso- 
lated dislocation in a single crystal of a purely random 
solutionJ^i2ii^ In the model Ni(Al) solution, the CRSS 
associated with an edge dislocation was found^^ to in- 
crease roughly linearly at a rate of about 30 MPa per 
atomic percent. This result agrees with the experimental 
data on the Ni(Al) hardness^ {H) from which the flow 
stress a can be deduced by application of the empirical 
linear rclation2Si21 H = 3<j (valid for metallic crystalline 



materials). The main question raised by such a result 
bears on the role of the screw dislocation segments. Ac- 
cording to a widespread belief drawn on the first order 
elastic dislocation theory (see text book as Ref. [28l ). the 
screw dislocation CRSS would be smaller than the edge 
one. By contrast, our present atomistic study shows that 
the screw CRSS has same order as the edge one for dif- 
ferent Al densities Cai, between 2 at. % and 12 at. %. 
Wc analyze our results in terms of the interaction be- 
tween the dislocation core and some isolated obstacles 
cither single Al placed at different positions around the 
glide plane or Al dimers with different positions, orienta- 
tions and bond lengths. The dominant component of the 
dislocation pinning is found to be a short range interac- 
tion of typically few Burgers vector between the obsta- 
cles and the distinct Schokley partials of the dissociated 
dislocation. This rather short range interaction is to be 
compared with the partial cores which spread over few 
lattice spacingsi ^^i^° Our systematic study of every dimcr 
configuration allows us to enhance the role of the chemi- 
cal interaction between nearest solutes which yields some 
pinning strengths that diverge from the simple superpo- 
sition of the strain fields due to each solute. 

In addition to provide some data about the strength of 
a single Ni(Al) crystal, our atomistic simulations allow us 
to challenge the different versions of the statistical the- 
ory for the SSHi^i^i^i^iiiSii^ The corresponding analytical 
models intend to provide an estimate of the CRSS from 
the elementary interaction between a single dislocation 
and an isolated obstacle. Most of the other effects on the 
dislocation pinning as the presence of the grain bound- 
aries, the dislocation forest and the thermal activation of 
the solute motion are assumed not to play an important 
role. This proves to match the conditions of our numeri- 
cal simulations which permits us to compare the theoret- 
ical predictions to the simulation data. A common point 
to the different versions of the SSH theory is to have been 
derived in the framework of the continuous line tension 
model where the dislocation is thought of as an elastic 
string anchored by a single type of obstacles. The split 
between the different versions of the theory stems from 
the various assumptions made on the critical string con- 
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figuration and how the dcpinning proceeds. For the sake 
of consistency, the string-obstacle interaction parameters 
are determined from our atomistic study and the statisti- 
cal models end results are compared to our direct compu- 
tation of the CRSS for the Ni(Al) single crystal. Among 
the theories proposed to compute the CRSS, we show 
that some of them quantitatively agree with our compu- 
tations for the concentrated solid solutions, i.e. 2 at. % 
< Cai < 12 at. %. The concordance of the atomistic 
simulations and the SSH continuous theory confirms the 
plausibility of a multi-scale approach to the plastic flow 
in the inhomogeneous media. 

The present paper is organized as follows: In Sec. II, 
the atomistic method to compute the CRSS for a screw 
dislocation is described and the results are compared to 
the edge ones. In Sec. Ill, the different analytical mod- 
els for SSH are presented and discussed in regard of our 
atomistic simulations. In Sec. IV, our study is resumed 
and our future works announced. 



II. MOLECULAR STATIC COMPUTATION OF 
THE SOLUTION STRENGTH 

By contrast to the phenomenological line tension ap- 
proach applied to the SSH,— the atomistic simulation 
allows to capture the main physical aspect of the dislo- 
cation core by integrating the nonlinear many-body in- 
teractions between the atoms displaced during the dislo- 
cation course. In our simulations the atomic interactions 
are modeled through the embedded atom method (EAM) 
the detail of which has been published elsewhere F^^i^i^i 
One must notice that two typo errors must be corrected 
before implementing the EAM potentials. For the Ni- 
Ni interaction see Ref. [H and for the Ni-Al interaction 
the coefficients gjqi and gAi must be exchanged in Ref. 
[T9I . According to private communications with other au- 
thors these corrections have been taken into account in 
other earlier works using the same method. The simula- 
tion cell (see FigdJ is oriented so as that the horizontal Z 
planes are the (111) of the face centered cubic (fee) lattice 
while the Y direction corresponds to the screw Burgers 
vector b — [110]ao/2 and the X direction is orthogonal 
to Z and Y and points at the dislocation motion. The 
simulation box size along the directions i = X,Y, Z is 
denoted by Li. The periodic boundary conditions are 
applied along X and Y while the external applied stress 
Ty2: is produced by imposing extra forces to the atoms 
in the upper and lower Z surfaces.— In order to form 
a screw dislocation between the two (111) central mid- 
planes, the displacement field of the elastic solution for 
a dissociated screw with Burgers vector b is applied to 
the atoms of the simulation box. In order to compensate 
the Burgers vector shift at the crossing of the boundaries 
along X, the corresponding periodic boundary conditions 
are tilted from b alongside Y. The ideal solid solution is 
then formed by substituting randomly Ni atoms with Al 
in the proportion fixed by the Al atomic density Cai- 




FIG. 1; (color online) View of the (111) glide plane and the 
dissociated screw dislocation in the Ni(Al) simulated solution 
at various concentrations: Cai — 2 at. % (Ihs), Cai = 6 at. 
% (center) and Cai = 10 at. % (rhs). The box size is 300 b 
along Y, 32 b along X and 34 b along Z (orthogonal to this 
paper sheet). The Shockley partials are colored in violet and 
the Al atoms in blue. The external shear stress has been fixed 
to one half of the CRSS. 



The solute distribution depends on the seed of the nu- 
merical random generator. For each Cai ranging from 2 
to 12 at. %, several distributions have been generated. 
The Molecular Statics (MS) simulations are performed to 
minimize the total enthalpy under a fixed applied shear 
stress. The external applied stress is incremented by 0.3 
MPa and for each increment the minimization procedure 
is repeated until it either converges to a required preci- 
sion or the dislocation starts to glide. To minimize the 
size effects in the simulations of the random solution, 
we chose Ly = 3006 which allows to neglect the inter- 
action between an obstacle and its periodic image and 
to achieve satisfactory statistics. The same method with 
same inter-atomic potentials have been employed in Ref. 
[25I for the edge dislocation. 

In Fig[T] we reported a snapshot of three typical sys- 
tems computed for an external shear stress smaller than 
the CRSS Tc. The atoms involved into the partials are 
recognized by their default configurations in first neigh- 
bor positions. The Al atoms that participate to the (111) 
planes that bound the glide plane have also been reported 
on these 3 pictures. The screw dislocation oscillates in 
the crystal in a way similar to the edge one as reported 
earlier.— It is worth noticing that such a wavy profile 
has been observed experimentally^ in some other type 
of solid solution as Cu(Al) and Cu(Si). The dislocation 
does not form large bow between well separated pinning 



points but rather conserves a wavy shape. Such con- 
figurations of the dislocation impinged on a random so- 
lution do not correspond to the one obtained within a 
phcnomenological elastic string model as the one used 
for instance in the work of Foreman and Makini^ The 
reason is threefold: (i) the solution is concentrated and 
thus the isolated point like obstacles arc rare, (ii) the 
interaction dislocation-obstacle is not a point-like force 
and (iii) the pinning strength of the Al obstacle is much 
weaker than what can be estimated through the elastic 
stress-strain field of a Volterra dislocation. As detailed 
further (see TabUl the pinning strength of an isolated Al 
on the screw dislocation is about one hundredth of the 
line tension meanwhile the elastic theory in a very first 
order version would predict one order larger. 

In our MS simulations when the dislocation moves 
freely, we let it glide for ten passages in the simulation 
box in order to probe possible stronger pinning configu- 
rations created by the relative displacement of one Burg- 
ers vector b at each passage between the plane above 
and below the glide plane. Since w 32 b, we assume 
that when the dislocation has glided over 10 x i^; we 
reach the maximum of the applied shear stress for the 
considered random distribution. Such a glide distance 
has the same order as the mean free path of a gliding 
dislocation through the forest dislocations in a metallic 
polycrystal where the dislocation density may reach 10^ 
cm~^. In most of our computations the dislocation does 
not encounter a new pinning configuration after 4 pas- 
sages in the simulation box. For each concentration, the 
calculations of the CRSS, Tc from various distributions 
have been reported in Fig[2] (a) where each open symbol 
corresponds to a different random distribution. The dis- 
persion on the measure of Tc is related to the finite size 
effect along Y. The choice of the suitable Ly results from 
a compromise between the computational load and the 
statistics. 

The order of magnitude of the solution strengthening 
is about 30 MPa per atomic percent of Al which is note- 
worthily similar to the edge one, obtained from Ref. [25l 
and reported in Fig|2] (b) for further comparison with 
the analytical models. Such a result contrasts seriously 
with the classical calculations based on the elastic theory 
of dislocationi^ However this can be fairly understood 
from the analysis of the interaction between the dislo- 
cation and a single obstacle at the atomistic level. We 
carried out the same type of MS simulations as for the 
random solid solution except that only one isolated Al 
atom is placed in a simulation cell of pure Ni. Then the 
external shear stress Tyz is incremented from zero to Tm at 
which the dislocation liberates. Since the simulation cell 
is periodic along Y, the obstacle and its periodic images 
form a regular array of obstacles separated by Ly. The 
balance between the Peach-Kohler force and the obstacle 
pinning leads to fm = (^m — Tp)bLy where Tp is the screw 
Peierls stress and b is the Burger vector of the whole 
dislocation. Because of the nonlinearity of the atomic 
interactions, we found that the pinning strength of an 
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FIG. 2: Variation of tlie critical resolved shear stress Tc 
for a screw dislocation (a) and for an edge dislocation (b) 
against the Al concentration Cm computed from the MS sim- 
ulations with different Al random distributions (symbols). 
The estimations made with the analytical models have also 
been reported: Fleischer- Friedel^i (dotted line), Friedel-Mott- 
Suzuki^ (full line), Mott-Nabarro-Labuscb^ (dot-dashed line) 
and Butt-Feltham^ (dashed line). 



isolated obstacle depends on its place above or below the 
glide plane as well as on which partial is concerned. Our 
results arc reported on the first lines of Tab|T] where the 
corresponding strength /„ has been normalized by the 
constant [il? to provide the standard pinning coefficient 
denoted by a. We choose [l = 74600 MPa, the shear mod- 
ulus for the Ni (111) planes [cn — cn -t-C44]/3. Although 
the normalization has been realized with the value of 
for pure Ni, one must bear in mind that this normaliza- 
tion is a conventional way to present the result. In Tab|T] 
the single obstacle denoted by (al) corresponds to an iso- 
lated Al placed in the (111) plane situated just above the 
glide plane while (bl) is for an Al which participates in 
the (111) plane just below the glide plane. The pinning 
strength of a single Al atom is found to have the same 
magnitude as for the edge dislocations^ which confirms 
our results about the same solid solution CRSS for both 
edge and screw dislocations. The strengths of some iso- 
lated Al have also been reported in TablH for the farther 
(111) planes. For the second nearest one, the pinning 



forces are referenced by (a2) and (b2) for the Al above 
and below the ghde plane, respectively. For the third 
nearest planes, the pinning forces are arranged in same 
order and referenced by (a3) and (b3). 

On the condition that the dislocation core distance to 
the obstacle remains larger than the core extent, it is pos- 
sible to analyze the interaction in terms of the continu- 
ous elastic theory. We intend applying such an approach 
for the obstacles situated in the third (111) plane from 
the glide plane. For that case, the potential landscape 
of the dissociated screw dislocation (dot-dashed line in 
Fig. O has the form l3z[l/{z^ + (x + df) - l/{z^ + x^)], 
i.e. the one for two opposite edge dislocations sep- 
arated from d, the dissociation distance i^'^'^^ Accord- 
ing to the elastic theory, the interaction pre-factor /3 is 
given by ^^j^(Ay) where v is the Ni Poisson coef- 
ficient. Ay is the atomic volume variation due to the 
Al impurity in Ni. If one neglects the interaction be- 
tween the obstacle and the farthest partial, the corre- 
sponding pinning coefficient for a single partial is thus 
a = i±^(Ay)6*y3/(87r62z2). The distance Z between 

the third (111) plane and the glide plane is Z = 56/-\/6. 
This corresponds to two and a half of the inter-plane dis- 
tance along the [111] direction. To provide an estimate 
of the pinning force we choose b* — h/^/Vl which is the 
edge component of a perfect Shockley partial. The vol- 
ume variation can be estimated according to a method 
described in the textbook^ AF = SvNiaQ^dao/dCAi, 
where the Vegard's law for Ni(Al) is used to express 
a^^dao/dCAi = 0.0763. Then one found a = 0.0015 
which is of same order though still larger than what has 
been found in our simulations (see lines (a3-b3) in Tab|I]). 
If instead of b* = b/ -\/l2, one chooses the effective Burg- 
ers vector computed through a Peierls-Nabarro method^S 
then a = 0.00083 which agrees better. A more quanti- 
tative study would require the account for the modulus 
misfit. This comparison allows to emphasize that as a 
consequence of the edge Burgers components of both par- 
tials the size effect dominates the modulus misfit effect 
often invoked in the analysis of the screw-impurity inter- 
action. Since the partials have some Burgers vectors that 
are not purely screw, a hydrostatic stress field is local- 
ized around the partial cores^ which renders the partials 
of a screw similar to the ones of an edge, at the atomic 
scale. Far away from the center of mass of the disloca- 
tion, since the edge components of the leading and the 
trailing screw partials are opposite, the hydrostatic stress 
field of both partials annihilate each other and the elastic 
theory prediction for a purely screw dislocation is recov- 
ered. In the first and second (111) planes, the shape of 
the interaction potential (full and dashed lines in Fig. [3l) 
is imposed by the nonlinear atomic interaction involved 
during the passage of the impurity through the core of 
the dislocation. Concerning the pinning strength on a 
single Al atom, we distinguish a common trend for screw 
and edge dislocations:^^ The anharmonicity enhances the 
pinning strength in the compressive regions in regard of 
the tensile ones. In our simulation cell, the compressive 
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FIG. 3: The internal potential energy for a dissociated screw 
dislocation against the dislocation obstacle distance. The ob- 
stacle is an isolated Al in the (111) planes above the glide 
plane: the first plane (full line), the second plane (dashed 
line) and the third one (dot-dashed). 



region of the edge dislocation is situated above the glide 
plane whereas for the screw 2 compressive regions can be 
distinguished, e.g. below the glide plane for the leading 
partial and above for the trailing partial. The obstacle 
labeled al (bl) in Tab. |T] visits the compressive region 
when it crosses the trailing (leading) partial. In agree- 
ment, the stronger pinning strength is found for the trail- 
ing (leading) partial. Furthermore the comparison of the 
pinning strengths in the separated compressive regions 
shows that the trailing partial is anchored more strongly 
than the leading one. The same trend can be noted for 
the tensile regions. Worthily the previous remarks are 
consistent with the results on the edge dislocationi^^ The 
latter trend finds some substantiations into the fact that 
the more stable position of the Al solute is inside the 
stacking fault ribbon, as can be seen in Fig. [3l The po- 
tential energy measured with respect to the isolated Al 
far away from the dislocation core should yield a diffu- 
sion current toward the dislocation stacking fault. This 
is a Suzuki-type effect which operates on both types of 
dislocation. This could lead to a classical dynamic strain 
aging of the fee substitutional alloys4S In our compu- 
tation the diffusion is frozen and such an effect is thus 
disregarded as it is indeed in the SSH statistical models 
discussed below. Another consequence of the absence of 
the solute diffusion is that there is no short range order- 
ing even for the Al concentrated solid solutions. It must 
be noticed that the interaction potential with the near- 
est obstacles (see Figl3]) is far from the simple shape as- 
sumed in some more phenomenological approaches^ and 
even from an elastic theory^i^ which neglects the dis- 
location core geometry. As it is required for the con- 
struction of a statistical model, some of the details of the 
atomic scale must be left off and one must retain only 
the essential quantities that are identified as important. 
According to the SSH theory, an obstacle can be charac- 
terized by an interaction range w and a pinning force /„, 
or equivalently a binding energy U. We thus attempt to 
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estimate such quantities and to that purpose we chose 
arbitrarily a rounding method. We distinguish only two 
different cases depending on the obstacle position with 
respect to the leading and the trailing partial. For the 
pinning strength a we choose the maximum of the inter- 
action force. The definition of the interaction range w 
requires a more tactful treatment since in principle the 
short range obstacle-core interaction is superposed to a 
long range elastic one and further the interaction poten- 
tial is not symmetric for an obstacle position ahead and 
behind the partials, mainly because of the presence of a 
stacking fault. Around the absolute force maximum, we 
measure w as the shortest distance to which the force van- 
ishes or falls to a local minimum. We choose to neglect 
the variations of the force over that limit. We concede 
that such a point can be discussed and must be kept in 
mind for further discussion on the analytical models. We 
emphasize that we have been primarily concerned with 
finding a manner to estimate one of the key parameters 
of the theory, w and that there is no well prescribed the- 
oretical way for that. We believe that the SSH theory 
should be developped to address specifically the case of a 
dissociated dislocation with an interaction potential sim- 
ilar to the one reported in Figl3] rather than the standard 
Gaussian-like potentials. 

In the present study, when the interaction potential 
shows a force maximum near the position of one of the 
partials we ascribe the force to this partial. This is a con- 
ventional way of arranging the numerous obstacle forces 
since the two partials are actually interacting through 
the stacking fault. Further when the obstacle is attrac- 
tive for the leading partial and repulsive for the trailing 
one, the force fields overlap and it may then be difficult 
to identify with accuracy the interaction range. When 
it has not been possible to separate different force max- 
imum, we reported the corresponding pinning strength 
as being exerted on the trailing one (see (b2-b3) lines in 
Tab|T|. In the fee symmetry, for the Al concentration 
we are concerned with, the density of the Al dimers has 
same order as the density of single Al, i.e. Cai ~ '^C'l; 
where n = 12 is the number of nearest fee neighbors. 
Above Cai — l/n the number of isolated Al vanishes in 
average. It is thus of some interest to study the pinning 
strength of the Al dimers that might be expected to play 
a role on the SSH because of the alloy ordering energy. 
Among the different configurations of pairs, we selected 
those with a distance between solute atoms correspond- 
ing to first, second and some of the third neighbors in 
the fee lattice. Either the screw dislocation interacts 
with preexisting Al dimers referred to as (c-h) in Tab[I] 
or the dislocation passage modifies the Al-Al bond cross- 
ing the glide plane (i-p) in TabUl For the non-crossing 
pairs, the (c-e) configurations correspond to the planar 
dimer situated above the glide plane whereas the (f-h) 
are below. The directions of the dimer bond before the 
dislocation passage are [Oil] (c) and (f), [101] (d) and 
(g), [110] (e) and (h), [Oil]. The dimers that cross the 
glide plane are oriented in the direction [721] (i), [211] 



TABLE L Summary of different pinning obstacles for both 
the leading (subscript I) and the trailing partials (subscript 
t) of the screw dislocation, their pinning force a (normalized 
by fib^) and their force range w. 



Nature 


text Ref. 


Qi and wi 


at and Wt 


single 


(al) 


0.0020/ 1.146 


0.0111/ 0.886 


1st planes 


(bl) 


0.0088/ 1.056 


0.0059/ 2.206 


single 


(a2) 


0.0036/ 1.416 


0.0032/ 1.696 


2nd planes 


(b2) 


- 


0.0053/ 1.526 


single 


(a3) 


- 


0.0011/ 2.86 


3rd planes 


(b3) 


- 


0.0009/ 2.856 


1st neighbor 


(c) 


0.0022/ 0.776 


0.0264/ 2.196 


non-crossing pair 


(d) 


0.0006/ 0.16 


0.0139/ 1.496 




(e) 


0.0033/ 1.036 


0.0168/ 1.146 




(f) 


0.0079/ 0.726 


0.0086/ 2.826 






0.0207/ 1.86 


0.0068/ 2.036 




(h) 


0.0153/ 1.056 


0.0106/ 2.336 


3th neighbor 


(i) 


0.0094/ 1.246 


0.0147/ 2.456 


crossing pair 


(j) 


0.0026/ 0.436 


0.0206/ 2.096 


2nd neighbor 


(k) 


0.0137/ 1.246 


0.0142/ 2.456 


crossing pair 


(1) 


0.0184/ 2.876 


0.0224/ 2.386 




(m) 


0.0096/ 1.096 


0.0146/ 1.076 


1st neighbor 


(n) 




0.0157/ 1.796 


crossing pair 


(o) 


0.0171/ 1.876 


0.0035/1.146 




(P) 


0.0218/ 2.776 


0.0161/ 1.376 



(j)_, [212] (k), [122] (1), [221] (m), [411] (n), [114] (o) and 
[141] (p). Our results on the pinning strength and the 
interaction range are reported on the corresponding lines 
in TablH Besides the variable geometry of the crossing 
pairs, it is worth noting that the pinning strength of a 
non-crossing dimer docs not correspond to the simple su- 
perposition of the strength of 2 isolated Al that would 
be placed at the reticular sites occupied by the dimer. 
The resulting strength can be even smaller than an iso- 
lated Al as seen from the comparison between the single 
(bl) and the first neighbor dimer (f). It must be also 
remarked that some of the dimer strengths can be larger 
than two isolated Al at the same position, e.g. (bl) and 
(g) in Tabdl The second neighbor Al dimers could be ex- 
pected to be more stable than others with regard to the 
LI2 ordering trend of the Ni(Al) alloy. Such an expecta- 
tion is not refiected by the dimer pinning strengths which 
are not significantly larger to break the second neigh- 
bor pairs although the interatomic potential between Ni 
and Al particles was adjusted to fit the corresponding or- 
der energyii^ The forces reported on lines (k-m) are not 
particularly stronger than those concerning the first and 
third neighbor dimers. Finally if one considers the whole 
set of the dimer configurations, the pinning strengths for 
a screw dislocation are not smaller than for the edge ones 
reported in Ref. [2^. 
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III. ANALYTICAL MODELS FOR COMPUTING 
THE SOLUTION CRSS 



The MS simulations allowed us to perform a direct 
computation of the CRSS at the atomic scale. In the 
past, different statistical theories have been proposed to 
evaluate the CRSS of solid solutions. The common as- 
sumption of the various analytical models is that the 
CRSS can be derived only from a unique interaction of 
the dislocation with an isolated impurity in the glide 
plane. Although these statistical models were primary 
devoted to the computation of the CRSS for an edge dis- 
location, there is no theoretical argument against the ap- 
plication of these models to a screw dislocation provided 
that one is able to quantify the model input parameters. 
In order to compare consistently the model predictions 
to our data on the SSH, the model parameters are deter- 
mined from the average of the single Al pinning configu- 
rations (lines (al) and (bl) in Tab|T|. We obtained a sin- 
gle valued obstacle strength = a/ifo^ where a = 0.007 
is the mean pinning coefhcient for the isolated Al situ- 
ated in the nearest (111) planes. We also computed the 
typical interaction range as an average over the different 
isolated obstacle configurations w = 1.36. Further, since 
in principle the models apply to an undissociated dislo- 
cation, the leading and the trailing partials are assumed 
to be tightly bound. In the different models, the end for- 
mula is usually presented for a theoretical square lattice 
which leads to an atomic area of s = 6^. In our system, 
this is to be changed for s = 6^-\/3/2,— i.e. the atomic 
area in the glide plane. 

Another quantity required before going into the details 
of the models is the dislocation stiffness. In our simula- 
tion cell, we computed the dislocation line energy as the 
difference between the whole atomic potential energy of 
a simulation cell with a dislocation and the one of a per- 
fect single crystal with same geometry and same atom 
number. It was then possible to obtain the energy per 
unit length El of the straight dislocation in pure Ni. One 
must emphasize that this computation implies the stack- 
ing fault energy which is not modeled in the classical es- 
timation of the dislocation energy, i.e.. El = fib^ /2. For 
the edge dislocation, we obtained Elc — 0-95El while for 
the screw Els = 0.61El- These two quantities depend on 
the simulation cell geometry (Lx,Lz) in agreement with 
the logarithmic law derived from elastic theory. Since 
in the Ni(Al) solid solution (see FiglT]) the dislocation 
length depends on the Al distribution, it is not possible 
to compute the dislocation line energy as done in a pure 
material where the dislocation is straight. However, one 
can expect reasonably that El should vary as the elastic 
shear modulus fi and the lattice parameter b, according 
to the Vegard's law which leads to a variation of roughly 
1 at. % of Cai- Since such a variation of the dislocation 
stiffness has not been considered in the SSH theory we 
shall not account further for such a dependency although 
for the concentrated solutions it would be worth account- 
ing for it. To obtain the dislocation line tension we apply 



a well known result from the isotropic elastic theory^ 
which tells us that for the screw Tg = i?Ls(l + f^)/(l — i^) 
while for the edge Fg = i?Le(l — 21^). We note that the 
line tension increase with the cell dimensions will impact 
the dislocation roughness through the obstacles distribu- 
tion as it is predicted by some of the following models. 



A. Fleischer-Friedel theory 

One of the classical models to evaluate the strength- 
ening of the solid solution is referred as to the Fleischer- 
Friedel (FF) model.— The model's main assumption is 
that the dislocation line forms some bows^ between iso- 
lated point-like obstacles randomly distributed. The bow 
shape assumption contrasts with what can be seen from 
our simulations in Fig[T] The critical angle associated 
with the impurity strength is = fm/2Ts. The Friedel 
length is Lp = ^2sT s/ cfm where c is the obstacle con- 
centration which corresponds to c = 4CAi to account 
for the different configurations, i.e. above and below the 
glide plane and for both partials (Tab[l|. The CRSS is 
then given by the balance between the Peach-Kohler force 
and a meanfield regular array of obstacle T^bLp = /„. 
The FF end formula reads: 



j-3/2 r 
by/2sTl' 



(1) 



The corresponding plot of against the Al density is 
reported in FigH] (a) and (b) for the screw and the edge 
dislocation, respectively. The FF theory is found to un- 
derestimate the CRSS of our simulations. Foreman and 
MakinSiii^ showed that the Eq[T] overestimates by about 
10 % the true CRSS for some perfect point-hke obstacle 
random distribution so that the discrepancy between our 
simulation and the FF model predictions could not be 
compensated by such a correction. It is compelling that 
the FF model is not suitable to compute the CRSS in the 
concentrated random Ni(Al) solution. One must be re- 
minded however that the FF model proved reliable^ for 
some stronger obstacles as precipitates and for smaller 
densities. It would be interesting to study by MS the 
lower densities but then the finite size sampling of the Al 
distribution limits the statistics. 



B. Mott-Nabarro-Labusch theory 

The Mott-Nabarro-Labusch (MNL) theory has been 
built from different contributions i^iii^ii^i^i^ The solute 
dislocation interaction is assumed to be of a finite range 
w without presuming of the attractive or the repulsive 
character of the interaction. We introduce the parameter 
V ~ 2'w which agrees with the definition of Nabarro for 
the interaction range. The dislocation is still idealized by 
a perfect elastic string and its configuration is assumed 
to be a quasi-straight line which interacts with several so- 
lute atoms. Alongside the dislocation core, in a ribbon of 
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length 2L and width 2u, the number of atomic sites in the 
planes contiguous to the glide plane is 8vL/s. The ribbon 
extent allows to account for the attractive and repulsive 
parts of the interaction. Inside this ribbon a counting of 
the obstacle gives in average 2n = 8vLc/s. The segment 
is sustained to a mean restoring force which depends on 
the stiffness of the string. Mott^ showed that this force 
could be written as 2Lf^x/L'v^Ts where x is the seg- 
ment mean position and L' is the mean distance between 
two obstacles situated along the segment 2L. The dis- 
tance L' is fixed by AL'vc/s = 1. Nabarro assumed^ that 
the characteristic length L could be identified as the dis- 
tance above which the string Green's function vanishes: 
L = (L'^I^s/^/2/™)2/3. Following Mott and Nabarro,-^ 
the counting of the obstacles situated in front and behind 
the segment 2L leads to a total number of interactions 
2n which the average force is ±/m/2. These interactions 
yield a maximum fluctuation of \/2nfm/2 which must 
be equated to the external Peach-Kohler force 2LbTc to 
obtain the CRSS: 



1/3 



(2) 



Here the obstacles above and below the glide plane have 
been accounted for. If one assumes that the two partials 
are bound, one must replace c with 2Cai- The plot of 
the corresponding CRSS has been reported in Fig[2] and 
it is found to underestimate the strengthening for the 
concentrated solutions. The same results hold for the 
edge dislocation. It is noteworthy that the MNL model 
is commonly thought to provide a good description for 
the high densities. In the case of our Ni(Al) solution, our 
comparison contrasts with such a belief. However, as for 
the FF model at smaller densities it is not excluded that 
the agreement could be recovered. An extended version 
of the MNL model has been proposed elsewhere^^ to ten- 
tatively account for the different types of obstacles, i.e. 
isolated and clusters with different pinning strength f„i 
and interaction range w. Although this model gave us 
some satisfactory results for both the edge and the screw 
dislocation further developments are required to be fully 
consistent on the theoretical treatment. 



C. Butt-Feltham theory 

In the Butt-Feltham (BF) theory^ i^^i^^ an undissoci- 
ated dislocation liberates from a pinning configuration 
by nucleation of a bulge. This critical bulge can be ap- 
proximated by a triangle shape of height W which cor- 
responds to the saddle point energy required to unzip 
the whole dislocation. According to Feltham W is esti- 
mated from the dislocation core radius extent and thus 
for a close-packed metals W « 3b. As other previous 
models in the BF theory it is proposed to relate the 
CRSS to the in-plane obstacle density. Along the dis- 
location line the inter-obstacle distance is roughly b/^/c. 



For a quantitative comparison we distinguish the mean 
distance in function of the different geometry for the edge 
and the screw dislocation: A = \f2bl\fc. and A = b/^/c^ 
respectively. The enthalpy required to form a bulge is : 
H = UL/\ + 2W'^T/L~'TbWL/2 where U is the binding 
energy and L is the bulge extent along the dislocation 
line. In the expression for i7, one recognizes the binding 
energy, the elastic cost for the bulge and the work of the 
Peach-Kohler force. The bulge curvature is determined 
by the external stress and thus for a certain W we deduce 
the corresponding bulge extent L ~ ^J8WT /rb. At the 
critical configuration, the enthalpy cancels which gives 
the end result: 



bW\' 



(3) 



We assume that the binding energy can be estimated 
roughly as ?7 = fmW and we account for the different 
obstacle configurations by c = ACai- The correspond- 
ing plot has been reported in Figsl2] (a) and [2] (b) for 
both edge and screw. Although the predicted CRSS is 
quantitatively of same order as the MS simulation data, 
the BF theory predicts a Tc in ^/Cai which underesti- 
mates the hardening rate of our simulations. It should 
be mentioned that Butt and Feltham also derived differ- 
ent power laws according to the variation of the bulge 
height W. With such an extension the BF theory pre- 
diction may agree better with the MS data. 



D. Priedel-Mott-Suzuki theory 

Another theory proposed by Friedel, Mott and Suzuki 
(FMS) in their text books^ considers the few atoms 
around an ideal undissociated dislocation line. This the- 
ory allows the dislocation to take locally much larger cur- 
vatures than in the FF model. The first effect of this 
pinning of the dislocation is to give the dislocation a 
"zigzag" shape. The interaction is characterized by a 
binding energy U which we approach by fmW valid for 
a linearized force. We briefly recall the model deriva- 
tion. The amplitude of the zigzag in the glide plane 
is denoted by W and its wavelength alongside the dis- 
location line is denoted by L. According to Friedel^ 
the number of solute contained in the glide plane area 
WL/2 is given by the inverse of the surface density, 
i.e. WL/2 = s/c. However a careful counting leads us 
to WL ~ s/c since the construction of a regular two- 
dimensional lattice where the obstacles are separated by 
L in the Y direction and W in the X direction leads 
us to WL ~ s/c which agrees with Suzuki's derivation^ 
The binding energy per unit of length is Ei = 2fmw/L. 
The line tension energy involved in the zigzag is given 
by E2 = r(V(W^2 + iV4) - L/2)/{L/2) « 2rwyL^ 
at first order in W/L. By minimizing the difference 
E2 — El one obtains the optimal value of the zigzag 
W = {fmWs/ATcY^^ that fixes the line shape and thus 
El = 2wWc/s and E2 = Wc'^W^/s^. The apphcation of 
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external shear stress deforms the dislocation. It has been 
assumed that the effect of the external shear stress is to 
unzip the line from the obstacles at the bottom of the 
zigzag, where the line tension exerts a maximum force 
in the glide direction. Considering that the maximum of 
the potential energy is a straight line {W = 0), bound to 
the obstacle situated at the top of the former zigzag, the 
energy difference between this maximum and the zigzag 
configuration is {—Ei/2) — {E2 — Ei). Multiplied by the 
wavelength L this gives the work that must be provided 
by the external shear stress to overcome the pinning bar- 
rier. The area that is comprised between the two config- 
urations of the line, zigzag and straight is WL/2, so the 
stress work is TcbWL/2 which must equal {Ei/2 — E2)L. 
The end result for Tc is: 



Noteworthily this expression for the CRSS does not de- 
pend on the line tension in contrast to the other models. 
To apply the theory to our system, the input parameters 
are estimated similarly and we assumed that c = 4Ca;. 
In Figl^Ja) and (b), it is remarkable that the Tc linear 
dependency in Cai is much closer from the simulations 
than a fractional power law. The FMS assumption of a 
zigzag shape is closer from the wavy profile of the dis- 
location seen in Fig[TJ The agreement between our data 
and the theory is better for the screw (Fig[2)^a)) than for 
the edge (FigOJb)) since for the latter the FMS model 
slightly overestimates the CRSS. With the force model 
parameters for a screw dislocation, the amplitude of the 
zigzag is = 0.436 at Cai = 2 at. % and W = 0.256 at 
Cai = 10 at. %. Besides the fact that such an amplitude 
is smaller than what can be depicted in Figll] the model 
predicts W to decrease with Cai ■ We studied the disloca- 
tion roughness and found the opposite trend. In view of 
the discrepancy on the edge CRSS, our procedure to com- 
pute the model parameters from our atomistic data can 
be discussed. To improve this transfer from the atomic 
scale to a continuous theory, we emphasize that some 
improvements of the model would be required to account 
for the dissociation of the dislocation core and the various 
types of obstacles. 

IV. SUMMARY AND PERSPECTIVES 

The aim of the present paper was to extend the study 
of CRSS for a solid solution Ni(Al) to a screw dislocation. 
We carried out two types of computation: (i) the CRSS 
of a random distribution of solutes at different concentra- 
tion Cai and (ii) the pinning strength of some different 
obstacles, isolated Al or dimers. We found that the CRSS 
and the pinning strengths of the screw dislocation are of 
the same order as the edge ones studied elsewhere.— Such 
a result could not be expected from a first order elastic 
theory which conventionally ascribes the screw pinning to 
a modulus misfit effect in contrast to the dominant size 



effect in the edge case. For the isolated Al impurities, 
this result has been clearly identified as a consequence of 
the dissociation of the dislocation core and the edge com- 
ponents of the Shockley partials cores. The Al clusters 
pinning forces have been shown not to derive from the 
linear superposition of the isolated obstacle strain field. 

Another issue of the present paper was to tentatively 
apply some analytical models to compute the CRSS. To 
that purpose, the elementary input parameters of the 
models were compiled from our atomistic data. The an- 
alytical models that we considered only account for the 
isolated foreign atoms situated in the planes that bound 
the glide plane. The Friedel-Mott- Suzuki theory provides 
us a better agreement for the Ni(Al) model solid solu- 
tion. The estimation of the interaction range between 
the dislocation and the obstacle is however to be ques- 
tioned because of the presence of the stacking fault and 
other nearest obstacles. Our work emphasizes that the 
SSH theory has to be developed to integrate more of the 
atomistic details of the dislocation solute interaction. We 
believe that the accuracy of the SSH theory in the fee 
solid solution requires to integrate the dissociation of the 
dislocation as well as the possibility for cluster forma- 
tion in the concentrated solution. The farther obstacles 
as those situated in the second and third (111) planes 
were found to yield a pinning strength that is not negli- 
gible in comparison to the solute atoms that bound the 
glide plane. We also tested other SSH theories involv- 
ing mixing laws as the Pythagorean one^ but they are 
essentially devoted to model the random distribution of 
point-like obstacles and finally did not allow us to obtain 
a better agreement. 

To predict quantitatively the CRSS in a variety of ma- 
terials, we believe that an extended use of the empiri- 
cal interatomic potentials is not feasible at the present 
state of our skill in developing such potentials. The re- 
cent progress of the first principles calculations for alloys 
and dislocations^Si^ allows us to expect that the pin- 
ning strength could be computed with a better accuracy. 
However it is difficult to imagine that the statistics could 
be also studied with such methods because of the as- 
sociated computational load. The effective interatomic 
potential can thus be employed to explore the frontier 
between the atomic details and the SSH statistics as we 
attempted in the present paper. In the near future, we 
shall propose to extend our work to the Al(Mg) solid 
solution. Our motivation rests on the fact that the inter- 
atomic potentials were proved physically reliable for the 
study of the alloy plasticit y^°i^'^ and the physical prop- 
erties of the Al(Mg) solid solution differ from the Ni(Al) 
ones, e.g. a larger size effect and a higher stacking fault 
energy. The temperature effect on the SSH and there- 
fore the dynamics of a dislocation in a disordered media 
is an eagerly difficult problem at the atomic scale because 
of the time limit of the molecular dynamics simulations. 
Moreover some puzzling problems of physics can be an- 
ticipated: (i) at the low temperature the inertial effect 
leads to a loss of strength)^ (ii) at higher temperature 
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the activation of the diffusion yields a dynamic aging^^ 
and (iii) for the intermediate temperature an odd ather- 
mal plateau on the yield stress has been clearly identified 
experimentally for various solutions and remains difficult 



to interpret through a simple theoryi^i^ We keep scruti- 
nizing the works about those topics which suffer from a 
lack of modern investigations. 
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